Face-Centered-Cubic B 8 o Metal: Density functional theory calculations 
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By means of ab initio calculations within the density functional theory, we have found that Bso fullerenes can 
condense to form stable face-centered-cubic (fee) solids. It is shown that when forming a crystal, Bso cages are 
geometrically distorted, the 1^ symmetry is lowered to T^, and four boron-boron chemical bonds are formed 
between every two nearest neighbor Bso cages. The cohesive energy of Bso fee solid is 0.23 eV/atom with 
respect to the isolated Bso fullerene. The calculated electronic structure reveals that the fee Bso solid is a metal. 
The predicted solid phase would constitute a new form of pure boron, and might have diverse implications. In 
addition, a simple electron counting rule is proposed, which could explain the stability of Bso fullerene and the 
recently predicted stable boron sheet. 

PACS numbers: 81.05.Tp, 61.50.Ah, 71.20.Tx, 74.10,+v 



Boron and carbon, the neighbors in the periodic table, both 
possess very rich physical and chemical properties. Car- 
bon forms the backbone of life, as there are millions or- 
ganic compounds that contain carbon in nature. Besides the 
graphite and diamond, carbon also possesses new isomor- 
phics, such as carbon nanotube and fullerenes [01 0], which 
have been extensively studied, and may be the fundamental 
materials for molecular electronics. Boron has a variety of 
complex isomorphic structures, such as a-rhombohedral B12, 
/^-rhombohedralBios, and tetragonal B50, etc (see e.g. OD for 
review). The boron solids are usually band insulators or semi- 
conductors S @1, but they could become metals or super- 
conductors under high pressure Q B In the past decade, 
boron sheet, cluster and boron nanotube (BNT) also gain wide 
attention , most of which were shown to be 

metallic, and they are expected to have broad applications in 
various circumstances. 

Recently, a new allotropic form of boron, Bgo, with a round, 
monoelemental, and hollow structure, has been predicted fl5ll . 
which is coined as boron fullerene. It consists of 80 boron 
atoms, and is very similar in shape to C60 fullerene except that 
an additional boron atom sits in the center of each hexagon 
(Fig. 1(a)). Bgo could be one of the most stable boron cages 
so far 11611 . which, if confirmed experimentally, should be 
the second example of a monoelemental buckyball after C60 
in nature. Our further analyses indicate that Bgo fullerene 
is probably formed by multi-center deficient-electron bonds, 
and satisfies a simple electron counting rule. This rule could 
also apply to explain the outstanding stability of the new class 
of boron sheets and nanotubes composed of triangles and 
hexagons predicted recently fl8n . On the other hand, re- 
call that soon after discovery of Ceo, people found that C60 
clusters can condense to form solid phases such as simple cu- 
bic (sc), face-centered-cubic (fee) lfl9H . and hexagonal-close- 
packed (hep) [20] crystals. As Bgo cluster has a geometri- 
cal structure similar to C60, and boron and carbon share some 
similarities in chemistry 112 ill , analogously, it would be natural 
to anticipate that the boron fullerenes could also condense to 
form solid phases. Here we show that, based on the ab initio 
calculations within the density functional theory (DFT) 12211 . 
Bgo clusters can indeed condense to form a stable fee solid 



that has a metallic electronic structure, in contrast to either 
the popular solid phases of pure boron that are usually insu- 
lating or semiconducting, or the fee and hep C60 that are band 
insulators [23]. In addition, we find that, unlike the case of 
C60 where the cage structure retains intact and is condensed 
by van der Walls force in solids JT9I Hi], when Bgo clusters 
condense to form the solid phase, the geometrical structure 
of Bgo is strikingly distorted [Fig. 1(d)], which is resulted 
from the formation of four boron-boron chemical bonds be- 
tween every two nearest neighbor Bgo clusters. The total en- 
ergy of Bgo fee solid is 0.23 eV/atom lower than the isolated 
Bgo fullerene, and 0.35 eV/atom higher than a-rhombohedral 
B12 solid. Therefore, if the fee Bgo metal is eventually con- 
firmed experimentally, it would compose a new form of pure 
boron in nature. 

Our calculations are mainly performed by means of the 
ABINIT package Il25n . This package is coded within the DFT 
framework based on pseudopotentials and plane waves, which 
relies on an efficient fast Fourier transform algorithm [26] for 
the conversion of wave functions between real and recipro- 
cal spaces, on the adaptation to a fixed potential of the band- 
by-band conjugate gradient method W27I1 . and on a potential- 
based conjugate-gradient algorithm for the determination of 
the self-consistent potential [Hsl E^l l30ll . Troullier-Martins 
norm conserving pseudopotentials [31] generated by fhi98PP 
code [32] are applied to mimic the electron-ion interaction, 
and the Perdew-Wang 92 1133m exchange-correlation potential 
within local density approximation (LDA) 113411 is used. The 
kinetic energy cutoff in the plane-wave basis is taken as 25 
Hartree, and the tolerance for absolute differences of the total 
energy is set as 10 -6 Hartree. 

The structure of the isolated Bgo cluster is optimized in a 
cubic supercell with a lattice parameter 25 A by the Broyden- 
Fletcher-Goldfarb-Shanno (BFGS) minimization l35ll with a 
convergence tolerance of 10 -5 Hartree/Bohr. As shown in 
Fig. 1(a), Bgo cluster, bearing 1^ symmetry, is similar in 
shape to C60, but there is an additional B atom centering in 
each hexagon, in agreement with the early calculation II 1 311 . 
Apparently, a pure two-center two-electron (2c2e) bonding 
scheme could not give a proper description of Bgo fullerene, 
as in which each boron atom has five or six nearest neighbors 
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FIG. 1: (Color online) (a) Schematic structure of an isolated Bso 
cluster; (b) is the view of green (dark) parts (a pair of neighbor 
hexagons) in (a) after a tt/6 rotation around the x axis; (c) the contour 
map of valence electron density in the plane of a hexagon in (b); (d) 
the distorted structure of Bso in the fee solid; (e) is the view of green 
(dark) parts in (d) after a tt/6 rotation around the x axis, showing 
the distortion of a pair of neighbor hexagons in (b); (f) the contour 
map of valence electron density in the plane of the rhombus, which 
is composed of the atoms marked by number 1-4 in (e). The red and 
blue colors in (c) and (f) represent high and low electron densities, 
respectively. (See the color bar in Fig. 2.) 

but it only possesses three valence electrons. We believe that 
the picture of the multi-center deficient-electron bonding ll36ll 
may apply to the present case. One may note that a hexagon 
with an additional boron atom sitting in the center can be 
viewed as a group of six triangles. A Bso fullerene consists 
of 20 such hexagons and 12 pentagon holes, or equivalently, 
120 triangles and 12 pentagon holes. Suppose that three boron 
atoms on the vertices of a triangle share a three-center two- 
electron (3c2e) bond. Then, a triangle would consume two 
electrons. As there are 120 triangles in Bso, 240 electrons are 
demanded for bondings. While the boron atoms on the ver- 
tices of pentagon holes have been counted in the triangles, no 
electron is consumed for pentagon holes. Thus, the total num- 
ber of demanded electrons is 240, which is just equal to the 
total number of electrons that Bso possesses. Based on the 
multi-center deficient-electron bonding picture, such a simple 
electron counting method may be generalized as follows: in a 
stable boron quasiplanar structure (sheet, nanotube, fullerene, 
etc.) composed by triangles and pentagon (hexagon) holes, if 
we assume that each triangle consumes two electrons, while 
pentagon (hexagon) holes consume no electrons, we would 
find that the number of total electrons demanded for bondings, 
that should be equal to the total number of valence electrons 
of the system, is nothing but twice the number of triangles. 



This electron counting rule also works for the newly predicted 
most stable boron sheet and nanotube (nlUE], which are com- 
posed by triangles and hexagon holes. In those cases, 8 boron 
atoms form a unit cell, which possess 12 triangles, so accord- 
ing to the above rule, 24 electrons in total are required for 
bondings, which just corresponds to the 24 valence electrons 
of 8 boron atoms, giving rise to a bonding balance. As illus- 
trated in Ref. 1 17], a pure hexagonal boron sheet is prone to 
accepting electrons, while a pure triangular boron sheet has a 
surplus of electrons, leading to that a mixture of these two 
phases with a proper proportion of triangles and hexagons 
would make the structure more stable. Our electron counting 
rule may be understood in a similar way. In a boron quasi- 
planar structure composed by triangles, a hexagon (pentagon) 
hole could be produced by removing the central atom from a 
group of six (five) triangles that share the central atom. As dis- 
cussed above, according to the multi-center deficient-electron 
bonding scheme, each triangle consumes two electrons, such 
a removal will reduce one boron atom that has three valence 
electrons, and therefore twelve (ten) electrons for bondings 
are cut down simultaneously. In this way, one may change 
the number of hexagon (pentagon) holes to realize the adjust- 
ment of the balance of the supplied and demanded electrons. 
Consequently, the electron counting rule manifests itself that 
a proper proportion of hexagon (pentagon) holes and triangles 
should be kept to give rise to the most stable boron quasiplanar 
structures. For instance, in Bso cluster, the proportion of pen- 
tagon holes and triangles is 1/10. Fig. 1(c) shows the valence 
electronic density of the hexagons in Bso cluster. One may 
see that high electron densities exist around the center boron 
atoms, showing a ring-like multi-center bonding pattern. As 
a hexagon with central atom [Fig. 1(c)] can be viewed as a 
combination of six triangles, the ring-like multi-center bond 
may be evolved from the combination of six 3c2e bonds. 

To seek for stable solid phases of B 80 , we have attempted 
sc and fee crystal structures, and taken two steps to find the 
optimized structures. In the first step, with the atomic posi- 
tions in Bso clusters fixed, the total energy of Bso solid with 
the presumed lattice, was calculated in every 0.1 i for the lat- 
tice constant (where the data with precision of 0.01 A are ob- 
tained by a spline interpolation from the calculated data with 
precision of 0.1 A). Without loss of generality, in the calcu- 
lations a properly fixed orientation of Bso molecules is kept, 
namely, the xy plane of Bso molecule is arranged to be along 
the (100) plane in the sc and fee structures, which makes the 
system bear a relatively high symmetry. In order to examine 
the effect of orientational disorder on the calculated results, 
we have rotated Bso clusters randomly to reduce the symmet- 
rical orientations, and recalculated the optimized lattice pa- 
rameters and the corresponding total energy. The results show 
that the total energy increase a lot (several eV per Bso cluster), 
suggesting that our presumed geometrical configuration is en- 
ergetically more favorable. We therefore get the total energy 
as a function of lattice parameters for both sc and fee struc- 
tures, where a single minimum in the curve of energy versus 
the lattice constant is found, indicating that a stable structure 
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may exist. In the second step, on the basis of the structures ob- 
tained in the first step, a full relaxation including the atomic 
positions, cell shape and volume was conducted by means of 
the BFGS minimization until the forces acting on atoms are 
less than a tolerance of 10 -5 Hartree/Bohr. As a result, after 
relaxation, the shape, lattice parameter, and the total energy 
are all varied, with the cohesive energies 4.9 eV and 18.2 eV 
per Bgo unit, i.e., 0.06 and 0.23 eV/atom, for the sc and fee 
Bgo solids, respectively, which are measured with respect to 
the isolated Bgo cluster. As the fee phase appears to be more 
stable in energy than the sc phase, we will focus on the fee 
phase in the following. 

As shown in Fig. 1(d), Bgo cage is dramatically distorted 
in fee solid, where the I/}, symmetry of an isolated Bgo is low- 
ered to Th symmetry for Bgo in fee phase. To check it further, 
we have relaxed a single distorted Bgo cage in a supercell, 
and found that it can indeed revert to the isolated Bgo cluster 
[Fig. 1(a)] with a spherical shape. Therefore, the distortion of 
Bso cages in a fee solid may be mainly owing to the formation 
of boron-boron chemical bonds between nearest neighbor Bgo 
cages. Fig. 1(b) and (e) show the configurations of a pair of 
hexagons along the z axis before and after distortion, respec- 
tively. The sharing arris of the two neighboring hexagons is 
elongated (from 1.678 A to 1.940 A) and broken, while two 
central atoms in the hexagons move outward, resulting in a 
planar rhombus structure composed by atoms denoted by the 
number 1-4 in Fig. 1(e), where the side length is 1.773 A 
and the angle is 57.87°. The rhombus is surrounded by four 
skew quadrilaterals, two triangles and two pentagons. Due 
to the symmetry, there are total six such pairs of hexagons 
along the x, y and z axes in one isolated Bgo cluster, which 
evolve into a deformed Bgo unit with six geometrically equiv- 
alent rhombuses. The boron atoms sitting in the centers of the 
rest eight hexagons in a Bgo unit all shift inward. Fig. 1(f) 
presents the valence electron density of the rhombus struc- 
ture. From the profile, we observe that the boron atoms la- 
beled by 3, 1,4 and 3, 2, 4 form two B-B-B bridge- type 3c2e 
bonds, which are very similar to the B-H-B bridge-type 3c2e 
bond in B2H6. The valence electron density is rather high 
between every two of the four boron atoms, and an elliptical 
high-density area presents in the center of the rhombus. It ap- 
pears that the two B-B-B bridge-type 3c2e bonds are coupled 
with each other, which eventually become a new complicated 
four-center bond. It will be illustrated later that the formation 
of this four-center rhombus is crucial to the bondings between 
the nearest neighbor Bgo units. 

As can be seen in Figs. 2(a) and (b), there are four boron- 
boron chemical bonds between every two nearest neighbor 
Bso units. The two green (dark) bonds are formed between 
two vertical rhombuses, while the two yellow (grey) bonds 
connect two parallel pentagons. The lengths of the green and 
yellow bonds are 1.734 A and 1.700 A, respectively, both 
falling into the typical boron-boron bond length scales in Bgo 
fullerene [15] and other boron materials. In a fee Bgo solid, 
every Bgo unit has 12 nearest neighbor Bgo units, and every 
two nearest neighbor Bgo units are bonded by four boron- 




FIG. 2: (Color online) (a) Schematic map of the bondings of the 
nearest neighboring five Bso clusters in (001) plane, the lengths of 
the green (dark) and yellow (grey) bonds are 1.734 A and 1.700 A, 
respectively; (b) a slightly rotational view of (a), where four boron- 
boron chemical bonds between the nearest neighbor Bso units can be 
seen clearly; (c) the contour map of valence electron density of the 
fee Bso solid (a) in the (001) plane through the center of Bso units, 
corresponding to the cross section of (a). The red and blue color 
represent high and low electron densities, respectively. The unit of 
color bar is electrons/Bohr 3 . 

boron chemical bonds, 48 bonds are therefore formed between 
every Bgo unit and its 12 neighbors. In this sense, the fee Bgo 
solid can be regarded as a three-dimension network of the dis- 
torted Bso units. Fig. 2(c) presents the valence electron den- 
sities of the fee Bgo solid in the (001) plane through the center 
of Bso units. The five big circles indicate the inner hollows of 
the five Bgo units in Fig. 2(a). It can be observed that there 
are high electron densities between neighboring cages, indi- 
cating stronger interactions between neighboring Bgo units. 
In particular, the densities corresponding to the green bonds 
are rather high, showing a a-like bond character. (The yellow 
bonds do not lay in this plane and are not shown directly in the 
density map.) Thus, we may conclude that Bgo clusters in the 
fee phase are connected by strong B-B a-like chemical bonds. 

In Fig. 2(c), it seems that small spherical-like hollow struc- 
tures are formed in the space between the neighbor Bgo units. 
By looking into it further, we find that it is indeed a B24 cage- 
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FIG. 3: (Color online) (a) The cage-like B 2 4 structure between the 
neighbor Bso units; (b) is the rotational view of (a); (c) is the dis- 
torted structure after relaxation in a supercell, where the rhombus 
breaks into two triangles. 



like structure composed of 6 rhombuses and 8 hexagons [Fig. 
3(a) and (b)]. This structure could be considered as a dis- 
torted truncated octahedron or cubo-octahedron, and interest- 
ingly, it is also a regular part of the crystal structure of metal- 
dodecaborites MB 12 (M is a metal) [37], in which a metal 
atom is located at the center of each boron cubo-octahedron. 
We have also relaxed the B24 cage in a supercell, and seen that 
it is distorted and becomes more spherical, and every rhombus 
breaks into two triangles [Fig. 3(c)]. The total energy of the 
relaxed B 2 4 cage is 0.57 eV/atom larger than the isolated B 80 
cluster. Thus, this B24 cubo-octahedral structure could exist 
only in a solid, and cannot appear in an isolated form. Similar 
to the situation of MB 12, the interstitial atoms may be inserted 
into these hollows to form new boron compounds. 

As listed in Table I, the total energy of fee Bgo solid is 0.23 
eV/atom lower than the isolated Bgo, 0.35 eV/atom higher 
than a-rhombohedral B12 solid. Note that the cohesive en- 
ergy of the fee C60 solid is about 1.6 eV per C60 L24], i.e., 
0.03 eV/atom. Furthermore, fee C60 solid and other molecu- 
lar solids J38, 39, 40] are condensed by van der Walls force, in 
which no chemical bonds are formed among the neighboring 
clusters. Our calculations on the inter-cluster bonds, valence 
electron densities, and the cohesive energies indicate that the 
condensed mechanism of fee Bgo solid is totally different from 
the above molecular solids. 

The energy bands and the density of states (DOS) for the fee 
Bso solid are calculated with the optimized lattice constants, 
as presented in Fig. 4. The energy bands are quite dispersive, 
and several bands spread across the Fermi level, as manifested 
in the left panel of Fig. 4, showing that the fee Bgo solid is a 
metal. As the most of the BNTs are also observed to have 
the metallic electronic structures, such a metallic property of 
the fee Bgo solid is conceivable. Unlike the fee C60 solid that 



is a band insulator with a direct energy gap of 1.5 eV 112411 . 



TABLE I: Total energies of the boron sheet, isolated Bso cluster, and 
Bso fee solid, with respect to a-rhomhedral boron solid, in eV/atom. 



boron sheet isolated B 8 o B 8 o fee solid 




Ref. [17] 
This work 



0.38 
0.46 



0.56 
0.58 



0.35 



FIG. 4: (Color online) Left panel: Energy bands of the fee Bso crys- 
tal. Right panel: Density of states (DOS) of electrons for the fee 
Bso crystal around the Fermi level, which is obtained by means of a 
smearing technique. The Fermi level is set to zero. 



the present result shows that the energy gap closes when Bgo 
clusters condense to form a fee solid. This is understandable, 
because there are strong interactions between the neighboring 
Bso clusters in a solid phase, which would inevitably enhance 
the overlap of the charge densities of electrons, thus broaden- 
ing the energy bands, and eventually leading to the vanishing 
of the energy gap. However, albeit the fee B 80 solid shows a 
metallic behavior, the electronic density at the Fermi surface 
is not so high, as shown in the right panel of Fig. 4, in which 
the DOS profile is depicted. The DOS has a minimum slightly 
above the Fermi level, giving rise to a low electronic density 
of states. Of particular interest is that the DOS profile of the 
fee Bso around the Fermi level look very similar to that of 
MgB2, in which the DOS of M gBg also exhibits a minimum 
just above the Fermi level iPH I^Il. MgB2 is a well-known 
two-gap superconductor with transition temperature T c ~ 39 
K 114311 . where it has been established that the energy bands at 
the Fermi level mainly derive from B orbitals, and is a typi- 
cal metal that is essentially attributed to metallic boron with 
covalent B-B and ionic B-Mg bonding |42|], whose supercon- 
ductivity is now well understood within the electron-phonon 
mechanism. This present analogy between the primary char- 
acters of electronic structures of fee Bgo and MgB2 gives a 
hint that superconductivity might exist in the fee Bgo crystal. 

In conclusion, in terms of the first-principles DFT calcula- 
tions, we have found that Bgo fullerenes can condense to form 
stable sc and fee solids. It is also uncovered that Bgo cages in 
fee solid phase are geometrically distorted, where the 1^ sym- 
metry is lowered to T^, and four boron-boron chemical bonds 
are formed between every two nearest neighbor Bgo cages. 
The total energy of Bgo fee solid is 0.23 eV/atom lower than 
the isolated Bgo fullerene. A simple electron counting rule 
is proposed, which could explain the stability of Bgo fullerene 
and the recently predicted most stable boron sheet. In compar- 
ison to the ordinary semiconducting boron crystals that can 
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become a metal or superconductor under high pressure, our 
calculated electronic structures, with some primary features 
quite similar to those of MgB 2 superconductor, show that the 
fee B 80 solid is a metal at ambient pressure, and may be a 
candidate of new superconductor. As a result, in spite of the 
popular phases of a-, (3- rhombohedral and tetragonal borons, 
the fee Bgo metal may be another novel form for pure boron 
in nature. 
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